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Abstract 

One-flavour QCD - a gauge theory with SU(3) colour gauge group and a fermion 
in the fundamental representation - is studied by Monte Carlo simulations. The mass 
spectrum of hadronic bound states is investigated in a volume with extensions of L ~ 
4.4 ro (— 2.2 fm) at two different lattice spacings: a ~ 0.37 rg (— 0.19 fm) and a ~ 
0.27 ro (— 0.13 fm). The lattice action is Symanzik tree-level-improved Wilson action 
for the gauge field and (unimproved) Wilson action for the fermion. 
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1 Introduction 



QCD with one flavour of quarks is an interesting tlieoretical laboratory to understand 
some aspects of the strong interaction dynamics, namely those not connected to sponta- 
neous chiral symmetry breaking and to the existence of light pseudo-Goldstone bosons. 
As a consequence of a quantum anomaly, the U(l) axial symmetry of the classical La- 
grangian is broken and in the limit of vanishing quark mass no massless Goldstone 
boson exists. 

An intriguing possibility at negative quark masses is the spontaneous breakdown of 
parity and charge conjugation symmetry - a phenomenon first conjectured by Dashen 
[T] in the three-flavour theory. This has to do with the possible negative sign of the 
fermion determinant at negative quark masses because under the assumption of the 
positiveness of the fermion determinant Vafa and Witten [2] proved the impossibility 
of this kind of spontaneous symmetry breaking. 

A dramatic consequence of the absence of (broken) chiral symmetry is the difficulty 
to find a unique definition of the point with zero quark mass in parameter space [5]. 
(For an excellent summary and discussion of this problem see [1].) 

Another line of recent theoretical developments is the relation between one-flavour 
{Nf = 1) QCD and supersymmetric Yang-Mills (SYM) theory with one supersymmetry 
charge {M = 1) [5j. This connection is the consequence of orientifold planar equiva- 
lence in the limit of large number of colours {Nc oo). This might imply approximate 
relations among hadron masses even at Nc = 3, for instance, the approximate degen- 
eracy of scalar and pseudoscalar bound states of quarks [B] reflecting the properties of 
the Veneziano-Yankielovicz low energy effective action of M = 1 SYM [7] in the mass 
spectrum of Nf = 1 QCD. For instance, the mass ratio of the lowest pseudoscalar 
meson to scalar meson is predicted, including 1/Nc corrections, to be {Nc — 2)/Nc [8]. 
Another prediction of orientifold equivalence is the size of the quark condensate in 
one-flavour QCD which has recently been compared with numerical simulation results 
in Ref. [9]. 

In the present paper we start to explore the mass spectrum of hadronic states in 
one-flavour QCD by numerical Monte Carlo simulations. This requires reasonably large 
physical volumes at small quark masses and high statistics - especially for determin- 
ing glueball masses and contributions of disconnected quark diagrams. We apply the 
Wilson lattice fermion action which has recently been shown by several collaborations 
[TUl [TH [n[ [T^ to be well suited for such an investigation. We start our exploratory 
studies here on 12'^ • 24 and 16^ • 32 lattices with lattice spacing a ~ 0.19 fm and 
a ~ 0.13 fm, respectively. This means that our present setup roughly corresponds to 
the earlier simulations of the qq-l-q Collaboration [lOj, but we hope to continue these 
investigations in the near future closer to the continuum limit as in Refs. [T H [T ^ [T ^ . 
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For setting the scale we use the Sommer parameter |14j tq which we set by definition 
to be ro = 0.5 fm. In other words, whenever we speak about "Ifm" we always mean 
"2ro" - having in mind that one-flavour QCD is a theory difl^erent from QCD realised 
in nature. 

Since the sign of the quark determinant is a sensitive issue, we carefully determine 
it and take it into account in determining the expectation values. In the present paper 
we choose the quark mass to be sufficiently far away from zero on the positive side, 
where the effect of the determinant sign is not very strong. In spite of this, as we 
shall see, we can investigate quite small quark masses down to rriq ~ 12MeV (that is 
nigro ~ 0.03), corresponding to a pion mass m-,^ ~ 270 MeV. 

Let us mention that keeping the quarks sufficiently heavy (choosing the hopping 
parameter k in the Wilson fermion action ([2]) below |) the problem of negative quark 
determinants can be avoided. (The thermodynamics of Nf = 1 QCD for heavy quarks 
have been investigated under this assumption in Ref. [15j.) Our aim is, however, to 
reach small quark masses and therefore we have to deal with the possibly negative sign 
of the quark determinant. 

For interpreting our results on the mass spectrum we find it useful to embed the 
Nf = 1 QCD theory in a partially quenched theory with more quark flavours. This 
embedding is particularly useful if the additional quenched valence quark flavours have 
the same mass as the dynamical sea quark because of the exact S\J{Nf) flavour sym- 
metry in the combined sea- and valence-sectors {Np denotes here the total number 
of quenched and unquenched flavours). In most cases we consider the natural choice 
Np = 3 which is closest to the situation realised in nature. We also work out some of 
the predictions of partially quenched chiral perturbation theory (PQChPT) and compare 
them to the numerical data. 

The plan of this paper is as follows: in the next section we define the lattice 
action and briefly discuss the updating algorithm. In Section [3] the partially quenched 
viewpoint is introduced and PQChPT is considered for it. Section [His devoted to the 
presentation of our numerical simulation data. The last section contains a discussion 
and summary. 

2 Lattice action and simulation algorithm 
2.1 Lattice action 

For the SU(3) Yang-Mills gauge field we apply, following Ref. [13| . the tree-level im- 
proved Symanzik (tlSym) action which is a generalisation of the Wilson plaquette gauge 
action. It belongs to a one-parameter family of actions obtained by renormalisation 
group considerations and in the Symanzik improvement scheme |16j . Those actions 
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also include, besides the usual (1 x 1) Wilson loop plaquette term, planar rectangular 
(1 X 2) Wilson loops: 

with the normalisation condition cq = 1— 8ci. For the tlSym action we have ci = —1/12 
The fermionic part of the lattice action is the simple (unimproved) Wilson action: 

^/ = E I - E [i^l+i.Uab,x,il + + - iM+f] I • (2) 

X i_i=l 

Here n is the hopping parameter related to the bare quark mass in lattice units amo 
by ^ 

— = amo + 4 . (3) 

The Wilson parameter removing the fermion doublers in the continuum limit is fixed 
in (ED-Q to r = 1. 



2.2 Simulation algorithm 

For preparing the sequences of gauge configurations a Polynomial Hybrid Monte Carlo 
(PHMC) updating algorithm was used, which is well-suited for theories with an odd 
number of fermion species. This algorithm is based on multi-step (actually two-step) 
polynomial approximations of the inverse fermion matrix with stochastic correction in 
the update chain as described in Ref. [18]. The starting point is the PHMC algorithm as 
introduced in Ref. [191 120j . The polynomial approximation scheme and the stochastic 
correction in the update chain are taken over from the two-step multi-boson algorithm 
of Ref. [2lJ . For details of the updating algorithm and for notations related to it see 
Ref. [E]. 

In order to speed up the updating even- odd preconditioning was used which pushes 
the small eigenvalues of the (squared Hermitean) fermion matrix to larger values. 

The eigenvalues of (5[C/]^ are assumed to be covered on typical gauge configurations by 
the approximation interval [e. A]. In exceptional cases some of the eigenvalues (typically 
just the smallest one) are outside this interval. In order to correct for this a correction 
factor C\U] is associated with such configurations. The exact value of this correction 
factor can be written as 

C[U] = { n [\'^^'""^ ^i(A.)^'2(A.)] y . (4) 

i 

Here the product runs over the eigenvalues of (5[C/]^, the polynomial ^1(2;) is an ap- 
proximation for 3;-i/(2"-s)^ P2{x) for [x^^^'^'^^^ Pi{x)\~'^ . The positive integer ub defines 



the determinant break-up which means that in the path integral the fermions are rep- 
resented by 

The part of the product in ^ where Aj is inside the interval [e, A] can be effectively 
replaced by a stochastic estimator and then 

N' 

c[u] = {]X[xy^''''^Pi{Xj)P2{Xj)]-^Y.^w{vU^-P\Q[un^^^ . (6) 

j n=l 

Here the Hj runs over the eigenvalues outside the interval [e, A], P'(x) is a sufficiently 
good approximation of [x^^^"^^^^ Pi{x)P2{x)]~^ , N' is the arbitrary number of stochas- 
tic estimators and the r]nS are Gaussian vectors in the subspace orthogonal to the 
eigenvectors corresponding to the eigenvalues Aj. In practice, one can choose the poly- 
nomial P2ix) to be such a good approximation that the stochastic part in ([6]) has no 
noticeable effect on the expectation values and therefore can completely be neglected. 
In this case the correction factor is simply given by 

C[U] = { n' [Af Pi{X,)P2{Xj)] . (7) 
j 

Besides the correction factor C[C/] , the sign (t[C/] of the fermion determinant det (^[t^] 
has also to be included in the reweighting of the configurations and then the expectation 
value of a quantity A is given by 

Jd[U]a[U]C[U]A[U] 
Jd[U]a[U]C[U] ■ 

This formula shows the dangerous sign problem which can arise due to the fluctuation 
of the determinant sign because in case of strong fluctuations of cr[U] both nominator 
and denominator on the right hand side may become small, spoiling the statistical 
accuracy. (Similarly, one can also loose statistics if the correction factors C[U] are 
much smaller than 1 on many configurations.) 

Typical values of the approximation interval and of the polynomial orders at the 
lightest quark mass simulated on 12'^ • 24 and 16^ • 32 lattices, respectively, are collected 
in Tabled! As in Ref. [18], the orders of the polynomials Pj, {j = 1,2) are denoted 
by Uj and those of Pj, {j = 1, 2) by fij, respectively. The simulations have been done 
with determinant break-up ub = 2. (The polynomials Pj are approximating {Pj)~2. 
For more details see |T8] and references therein.) 

The last four columns of Table [T] show the values of the deviation norm 5 which is 
minimised for a given polynomial order n in the least-square approximation scheme we 
are using. Generically 6 is defined as 

f L'dxw{x)[f{x)-P^{x)f Y2 
^ dxw{x)f{xY ^ 



Table 1: Algorithmic parameters in the runs with lightest quark mass on 
12^ ■ 24 (first line) and 16^ ■ 32 (second line) lattice, respectively. For 
notations see the text and also Ref. [18]. 



e 


A 


ni 


ni 


^2 


^2 


Si 


Si 


S2 


S2 


3.25 ■ 10-6 


2.6 


350 


550 


1400 


1600 


4.9 ■ 10-^ 


6.7-10-^ 


9.9 ■ 10-^ 


8.8 ■ 10-^ 


1.2 ■ 10-^ 


2.4 


250 


370 


1000 


1150 


5.4 ■ 10-^ 


8.2 ■ 10-^ 


4.8 ■ 10-^ 


3.1 ■ 10-^ 



Here /(x) is the function to be approximated and w{x) is a positive weight function 
actually chosen in our case to be wi{x) = W2{x) = x^/(^"s) g^^(^ wi{x) = W2{x) = 1, 
respectively. The values of 6i in Table [1] are such that the average acceptance rate of 
the stochastic correction at the end of trajectory sequences is between 80 — 90%. The 
other 6 values are small enough to ensure practically infinite precision of the expectation 
values. For more details on the algorithmic setup in our runs see also Section 21 

3 Partially quenched viewpoint 

Because the classical U(1)a axial symmetry is anomalous, the single-flavour QCD the- 
ory does not have a continuous chiral symmetry apart from the U(l) quark number 
symmetry. Consequently it does not have spontaneous chiral symmetry breaking and 
hence no (pseudo-) Goldstone bosons and no easy definition of the quark mass [3]. In 
the lattice regularisation it is, however, possible to enhance the symmetry artificially 
by adding extra valence quarks which are quenched, that is, are not taken into account 
in the Boltzmann- weight of the gauge configurations by their fermion determinants. In 
principle, one might consider any number of quenched valence quarks with any mass 
values but, to remain close to QCD realised in nature, the most natural choice is to 
take two equal-mass valence quarks and to call them u and d quarks. The original 
dynamical quark can then be called s quark where "s" may stand for sea or strange. 
The theory with dynamical s quark and quenched u and d quarks is partially quenched. 
(Observe that this partially quenching is somewhat unconventional, since some of the 
valence quarks are quenched but taken degenerate with the sea quark.) 

Using this terminology, for instance, the pseudoscalar bound state of s and s can be 
called r]s. The corresponding scalar state is then ag- The lowest baryon state consisting 
of s quarks, which has to have spin | because of the Pauli principle, can be named 0~ 
or e.g. As etc. 

A theoretical description of partially quenched QCD can be obtained through the 
introduction of ghost quarks [H]. For each (quenched) valence quark a corresponding 
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bosonic ghost quark is added to the model. The functional integral over the ghost 
quark fields then cancels the fermion determinant of the valence quarks and only the 
sea quark determinant remains in the measure. In our case there are 2 flavours of 
valence quarks and ghost quarks, each, with equal masses my, and a single flavour of 
sea quarks with mass ms- 

A particularly interesting point of the partially quenched theory is the one where 
all the three quark masses are equal. In this point there is an exact SU(3) vector-like 
flavour symmetry in the valence + sea quark sector, and the hadronic bound states 
appear in exactly degenerate SU(3)-symmetric multiplets. For instance, there is a 
degenerate octet of pseudoscalar mesons - the "pions" (vr'^, a = 1, . . . , 8) satisfying an 
SU(3)-symmetric PCAC relation. With the help of the divergence of the axialvector 
current and pseudoscalar density one can define, as usual, the bare PCAC 
quark mass ampcAc in lattice units: 

Here the indices + and — refer to the "charged" components corresponding to Xa ± iX), 
(with Xa^b some off-diagonal Gell-Mann matrices) and d* denotes the backward lattice 
derivative. Due to the exact SU(3)-symmetry, the renormalised quark mass correspond- 
ing to mpcAc can be defined by an SU(3)-symmetric multiplicative renormalisation: 



By tuning the bare quark mass on the lattice suitably, the masses of the "pions" can 
be made to vanish, as the numerical results indicate, and the renormalised quark mass 
vanishes, too. At this point the partially quenched theory has a graded S\J{Nf\Nv)l ® 
SU(A''ir|A^y)/j symmetry, which is broken spontaneously to a "flavour" SU(A''i7'| A^y). 
(Here Ny is the number of additional valence quark flavours and Np = Ny + Nj = 
Ny + 1.) In our case, with Ny = 2 flavours of valence quarks, the symmetry is thus 
SU(3|2). The "pions" are the Goldstone bosons of the broken SU(3) subgroup. 

Adding generic quark masses my and ms, the symmetry group is explicitly broken 
down to SU(2|2). In the special case my = ms, considered here, the symmetry is still 
SU(3|2), and its subgroup SU(3) is the flavour symmetry mentioned above. 

The "pions" are, of course, not physical particles in the spectrum of Nf = 1 QCD. 
Nevertheless, their properties such as masses and decay constants are well defined 
quantities which can be computed on the lattice. The same is true of the PCAC quark 
mass ^T^pcAC' which is therefore a potential candidate for a definition of a quark mass 
of this theory. 

The relation between the pion masses and the quark masses can be considered in 
partially quenched chiral perturbation theory [25\ [26], including effects of the lattice 
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spacing a [271 EEl EHI EOj EI]. The pseudo-Goldstone fields are parameterized by a 
graded matrix 

U{x) = exp (jrH^)^ (12) 

in the supergroup SU(3|2). (Here the normahzation of Fq is such that its phenomeno- 
logical value is ~ 86MeV.) The commuting elements of the graded matrix $ represent 
the pseudo-Goldstone bosons made from a quark and an anti-quark with equal statis- 
tics, and the anticommuting elements of $ represent pseudo-Goldstone fermions which 
are built from one fermionic quark and one bosonic quark. The supertrace of $ has to 
vanish, which can be implemented by a suitable choice of generators |32j . 

We have calculated the masses of pseudo-Goldstone bosons in next-to-leading order 
of partially quenched chiral perturbation theory along the lines of Ref. [32], including 
0(a) lattice effects [2^. The quark masses enter the expressions in the combinations 

Xy = 2^0 my, xs = '^Boms, (13) 

with the usual low-energy constant Bq, and the lattice spacing occurs as 

p = 2Woa, (14) 

where Wq is another, lattice-specific, low-energy constant. For the pion masses we 
obtain 



2 _ 2 , , XV + P 

myy = m^ = XV + P+ .r, 2t?2 



(2xy - + /?) In ( '^^J^p2 ) +XV -XS 



+-2 [(2L8 - Lr>)xv + (2^6 - U)xvXs 

+ {2W8 + Wq-W5-W^- L^)pxv + - U)PXS] , (15) 

where the usual low-energy parameters Lj appear, together with addtional ones (Wi) 
describing lattice artifacts. 

The mixed mesons, whose masses mys we have also calculated, become degenerate 
with the pions in the special case my = ms- In this case the expression reduces to 



2 ^ , ix + p)\ ( X + P \ 



+--2 [(2L8 - Ls + 2L6 - U)x^ 

^0 

+{2Wf, + 2W^ -W^-W^-L^- U)xp\ . (16) 



To leading order the PCAC quark mass obeys 2Bq m^^^^ = X + P-, and we recognize 
the Gell-Mann-Oakes-Renner relation 



=2Bom^cAc + NLO. (17) 
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Including terms in next-to-leading (NLO) order, we can express in terms of mpQ^^ 
as 

2 

2 _ , XpCAC 1 Xpcac 



+ J2 [(2^8 - ^5 + 2L6 - L4)Xpcac 







+{W^ + WQ-W^-Wi- 2L8 + L5- 2Lq + L4)xpcacP] , (18) 



where we define 



Xpcac = 2-Bo m^cAc • (19) 

As a remark, in the case niv = ms the masses can alternatively be obtained from 
the partially quenched theory with symmetry SU(2|1) by considering mixed pions made 
from a valence quark and a degenerate sea quark. Indeed, calculating the masses in 
this model reproduces (fT6|) . 

The rjs can be included in the analysis by relaxing the constraint of a vanishing 
supertrace [SSI [32], and associating it with the field 



4'o(x) = sTr$(x). (20) 

The effective Lagrangian then contains additional terms depending on $0- 

AC = ad^^od^^o + ml<^>l + 0{^l), (21) 

where a and m$ are free parameters in this context. We content ourselves with dis- 
playing only the leading order expression for the mass of the r]s , which reads 



ml = < + (22) 



Our numerical results for m^^ allow to determine a and m^. 



4 Numerical simulations 

After some preparatory search in the parameter space we concentrated our runs on 
the 12^ • 24 lattice to /5 = 3.8 and those on 16^ • 32 to /5 = 4.0. The parameter 
values, the number of analysed configurations, the average plaquette, its integrated 
autocorrelation and the value of the Sommer scale parameter in lattice units rg/a are 
summarised in TableO As one can see, taking the values of ro/a at highest k's (smallest 
quark masses), the extensions of the 12^ and 16^ lattices are L = 4.46 = 2.23 fm 
and L = 4.29 ro = 2.14 fm, respectively. Since we fix tq = 0.5 fm by definition, these 
correspond to lattice spacings a = 0.186 fm and a = 0.134 fm, respectively. 

In the update-chain by the PHMC algorithm with stochastic correction [18] a se- 
quence of PHMC trajectories is followed by a Metropolis accept-reject step with a 
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higher precision polynomial. The total length of the trajectory sequence in the runs 
in Table [5] was between 1.5 and 1.8. The sequences consisted out of 3-6 individual tra- 
jectories. The precision of the first step of polynomial approximations was tuned such 
that the acceptance of the PHMC trajectories was about 0.80-0.85. The total length 
of the trajectory sequence was chosen such that the acceptance of the Metropolis test 
was again 0.80-0.85. This ensured a relatively high total acceptance of 0.64-0.72. Dur- 
ing the runs we tried to optimise the parameters of PHMC. The different values of 
the integrated autocorrelation times for the average plaquette in Table [2] are, in fact, 
mainly due to increasingly better optimisations and not so much to the dependence on 
run parameters. 

The second step approximations were more than good enough to ensure that the 
expectation values were completely unaffected by the remaining small imprecision. 
(See, for instance, the small relative deviations in Tabled!) This has also been explicitly 
checked by performing a final stochastic correction on a large sample of configurations 
with polynomials P' of order 2500 in the stochastic part of the right hand side of ([6]). 

For the calculation of the expectation values the reweighting procedure according 
to ([8|) has to be carried out. For this, besides the correction factor C[U] from d?]), also 
the sign of the fermion determinant cr[U] is needed. This we calculated by the spectral 
flow method [22]. For the K-dependent computation of the low-lying eigenvalues of the 
hermitean fermion matrix Q[U] we followed Ref. [23]. 

It turned out that the effect of the correction factors cr[C/]C[f7] is in most cases 
negligible. For instance, in run b of Table [2] the average value of (t[C/]C[C/] in the 
denominator is 0.9982. In run c it is 0.9842. In run b there are 34 configurations out 
of 3403 where some eigenvalue is outside the approximation interval [e. A] and out of 
them there is a single one with negative fermion determinant. In run c there are 167 
from 2884 outside [e. A] and out of them there are 26 with negative correction factor 
due to 0" = — 1. 

Since the sign of the fermion determinant was not determined on every configura- 
tion, the question arises whether perhaps some negative signs were missed. This is very 
improbable because we determined the sign also on the neighbouring configurations in 
addition to those with small eigenvalues and out of the remaining configurations we 
have chosen 100 randomly for sign determination. None of these additional configura- 
tions turned out to have a negative determinant. 

In the average plaquette and rg/a the effect of the correction factors is completely 
negligible. For instance, in runs b and c the correction has an effect in the average value 
of ro/a only in the fifth digit - whereas the statistical error is in the third digit. In all 
other runs besides b and c every eigenvalue is inside the approximation interval [e. A] and 
therefore, according to ([7]), the correction factor is equal to 1 on every configuration. 
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Table 2: Summary of the runs: 12^ ■ 24 and IQ^ ■ 32 lattices have lowercase 
and uppercase labels, respectively. The number of gauge configurations, 
which were saved after every trajectory sequence, is Ai'conf- The average 
plaquette value, its autocorrelation in number of trajectory sequences Xpiaq 
and the value of r^/a are also given. 



label 






^conf 


plaquette 


^plaq 


ro/a 


a 


3.80 


0.1700 


5424 


0.546041(66) 


12.5 


2.66(4) 


h 


3.80 


0.1705 


3403 


0.546881(46) 


4.6 


2.67(5) 


c 


3.80 


0.1710 


2884 


0.547840(67) 


7.6 


2.69(5) 


A 


4.00 


0.1600 


1201 


0.581427(36) 


4.3 


3.56(5) 


B 


4.00 


0.1610 


1035 


0.582273(36) 


4.1 


3.61(5) 


C 


4.00 


0.1615 


1005 


0.582781(32) 


3.3 


3.73(5) 



4.1 Results for hadron masses 

Starting with the mesonic states, we consider the simplest interpolating operators in 
the pseudoscalar and scalar sectors: 

0+ : P(x) = i^ixhni^ix) , (23) 
0" : S{x) = i:{x)ip{x) . (24) 

We denote with ijg and ag the corresponding hadron states at the lowest end of the 
energy spectrum (the usual notation is used for the respective quantum numbers). 
Corresponding states in the QCD spectrum with the same quantum numbers are the 
r]'(958) and /o(600) (or a). (Note, however, that the states in QCD are linear combi- 
nations of uu, dd and ss components - in contrast to the states in Nf = 1 QCD which 
are built out of a single quark flavour.) 

In the case of the pseudoscalar mesons, invariance under the flavour group plays 
a special role when comparing with QCD states because of the U(l) axial anomaly. 
(This is not the case for baryons, see the following.) 

Analogously to flavour singlet mesons in QCD, the correlators of the above inter- 
polating operators contain disconnected diagrams. These were computed by applying 
stochastic estimator techniques (SET), and in particular the variant of [33] with Z2 
noise and spin dilution. The method was already applied to the case of SYM [3l] (as 
mentioned in the introduction, SYM shares many similarities with Nf = 1 QCD). In 



11 



Table 3: Results for light hadron masses in Nf = 1 QCD. 



run 


arrir, 

'Is 


(-* s 


amQ++ 


am^^ 


a 


0.462(13) 


0.660(39) 


0.777(11) 


1.215(20) 


I. 




U.4(J3(llj 


0.629(29) 


0.685(10) 


1.116(38) 


c 


0.398(28) 


0.584(55) 


0.842(16) 


1.204(57) 


A 

A 


n /I /'I 
U.400(^i 1 ) 


U.oU i{ol) 


i.(}oo[ ly ) 


i.UUo(^i5 ) 


B 


0.380(18) 


0.554(52) 


1.032(66) 


0.960(15) 


C 


0.316(22) 


0.613(67) 


0.980(97) 


0.876(26) 



order to optimize the computational load, taking also autocorrelations into account, 
every fifth configuration was typically analysed, with 20 stochastic estimates each. 

Spin states can be also build by purely gluonic operators. These are a well known 
object of investigation in lattice QCD were they should describe the glueballs. Due 
to the expected signal-noise ratio of their purely gluonic correlation they belong to 
the most notorious particles to measure. In particular the 0"*"^ glueball has the same 
quantum numbers as the as meson. As a consequence, these two states can also mix 
with each other but in this first investigation we neglect the mixing and consider only 
diagonal correlators for both states. 

We used the single spatial plaquette to obtain the mass of the 0"^"'" ground state. 
To increase the overlap of the operator with this state we used APE smearing and 
also performed variational methods to obtain optimal glueball operators from linear 
combinations of the basic operators. 

We now come to the baryon sector. The simplest baryonic interpolating field which 
can be built out of one quark flavour is 

Ai(x) = €abc['4'aix)'^C-fi1pbix)]lpcix) ■ (25) 

The above operator also contains a spin 1/2 component implying that the spin 3/2 
component, on which we focus, must be projected out from the spinorial correlator 

G,iit) = ^{Ajix,t)Ai{0)) . (26) 

X 

We follow [35j and consider the spin-projected correlator 

Gs/2{t) = ^Tr [Gji{t)j,j, + Gu{t)] . (27) 
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The low lying hadron state contributing to the above correlator is expected to have 
positive parity (|^). This corresponds to the A(1232)++ of QCD if our dynamical 
fermion is interpreted as an u quark. If the dynamical fermion is taken to be the s quark 
then this would be the 0,~ baryon. (However, spin and parity of the corresponding 
particle have not been yet measured, so the identification of this state with the Q~ 
baryon is still uncertain |36|). In correspondence to rjg and ag, in what follows we call 
this state A^. (Here one can interpret the index s as referring to the "sea" quark.) 

It should be noted at this point that the above QCD states are not flavour singlets 
in Nf = 3 QCD (and in the one flavor partially quenched theory). We recall here that 
interpolating fields corresponding to flavour singlet baryon states cannot be build in 
QCD if only quark fields are considered as ingredients. 

The results of the hadron masses are reported in Table [3] and, as a function of the 
bare PC AC quark mass ttt-pcao in Fig- E In the figure the masses are multiplied by the 
Sommer scale parameter rg, therefore one can put the results for both lattice spacings in 
a single plot and check their scaling. (The expected small change of the multiplicative 
renormalisation factor of mpcAc between /3 = 3.8 and /3 = 4.0 is neglected here.) 

Only in the case of run c the measurement correction has a sizeable effect on the 
mass estimates. In this case configurations with negative determinant where singled 
out: the sign of the determinant has the effect of pushing the masses up by 7 — 10 %. 

The errors on the glueball mass are rather large - especially on the 16^ • 32 lattice 
at /3 = 4.0 - therefore they are not shown in the figure. Obviously, our statistics is not 
sufficient for this purpose. In general a larger number of configurations would improve 
the determinations in the glueball sector. Since the computational load is in this case 
negligible, for future runs we plan a more frequent storage of the gauge configuration. 

4.1.1 Valence analysis 

The connected contribution to the meson correlators can be interpreted as a non sin- 
glet meson made up of valence quarks in the partially quenched picture, see Sec. El 
The pseudoscalar channel corresponds in particular to the "valence" pion. Since the 
computation of the connected diagrams is less demanding, we could afford the analysis 
of the complete set of configurations. 

In the baryon sector, one can define a "valence" nucleon, with the usual projector 
operator 

N{x) = eabciMxfCi^iixMci^) , (28) 

where ip' can be interpreted as the field of the valence quark. 

The results concerning valence hadron masses are reported in Table [H and Fig. [2j 
In addition, the bare PCAC quark mass according to the definition in (jlOp and the 
bare pion decay constant in lattice units af.,^ are also included, and its renormalised 
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Figure 1: The mass of the lightest physical particles in one-flavour QCD 
as a function of the PC AC quark mass. The masses are multiplied by the 
scale parameter Tq in order to obtain dimensionless quantities. Open and 
full symbols refer to f3 = 3.8 and (3 = 4.0, respectively. 



counterpart f^ are defined as 

af^ = (amO-i(0|.4+^o,/.=ok"(^^= 0)) , = ^a/. (29) 

where A^^ is the axialvector current as in (llOp and vr" (p = 0) is a pion state with zero 
momentum. (The normalisation of is such that in nature we have f^ ~ 130 MeV.) 
The value of a/^ on the lattice is obtained by the method described in |37j . In Fig. [2] 
the masses are multiplied by the scale parameter tq in order to obtain dimensionless 
variables. 



4.1.2 Chiral Perturbation Theory fits 

The properties of the valence pion (pion mass 771,^ and decay constant /^) can be 
analysed in partially quenched ChPT. We fit a'^m'^ and o/tt simultaneously as a function 
of 

ampcAc including the data at both values of /5. There are not enough data in order to 
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Table 4: The PC AC quark mass mpcAC; the pion mass and decay 
constant fj^, and the nucleon mass mj\f in lattice units. 



run 










a 


0.02771(45) 


0.3908(24) 


0.1838(11) 


1.0439(54) 


b 


0.01951(39) 


0.3292(25) 


0.1730(15) 


0.956(27) 


a 


0.0108(12) 


0.253(10) 


0.156(10) 


1.011(51) 


A 


0.04290(36) 


0.4132(21) 


0.1449(9) 


0.9018(44) 


B 


0.02561(31) 


0.3199(22) 


0.1289(10) 


0.7978(53) 


C 


0.01700(30) 


0.2635(24) 


0.1188(12) 


0.734(10) 



account for the lattice artifacts. Therefore the fit is done with the continuum formulae 

2 , XpCAC T XpCAC /tt 1 XpCAC , XpCAC /on\ 

-. = XPCAc + ^^ln^, ;^ = 1-^^1-^' (30) 

with the low-energy constants 

A3 = 47rFoexp{647r2(L4 + L5-2L6-2L8)}, 

A4 = 47rFoexp{647r^(L4 + L5)}. (31) 

The changes of the renormalisation constants Z^, Zp between the two (5 values are 
neglected. The results are displayed in Figs. [3] and [H 

Owing to the fact that the number of degrees of freedom in the fit is small, the 
uncertainty of the fit parameters is relatively large. The determination of the universal 
low-energy scales A3/F0 and A4/F0 can be improved by considering the ratios [iOl [TU] 

7^, (32) 

in which some of the coefficients cancel. We consider the data on the larger lattice at 
fi = 4.0 and take the quantities at k = 0.1615 as reference. The fit yields 

^ = 10.0 ±2.6, (33) 

^ = 31.5 ±14.3, (34) 

which is compatible with the phenomenological values from ordinary QCD [38]. 

In order to estimate the parameters a and m$, related to the mass of the rjs (see 
Sec. [3]), we made a fit of and m^^ at /? = 4.0 in leading-order ChPT. The result is 

a = -0.03(19) , am$ = 0.18(8) , (35) 
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N, (3=3.8 ^ 
7C, p=3.8 
N, p=4.0 
71, p=4.0 
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0.05 0.1 
''O'^PCAC 

Figure 2: The mass of the valence pion and nucleon as a function of the 
bare PCAC quark mass. Open and full symbols refer to (3 = 3.8 and 
(3 = 4.0, respectively. 



indicating the vanishing of a. Fixing a = in the fit yields 

am$ = 0.19(2) or rom$ = 0.72(10) , 



(36) 



where the value of ro/a extrapolated to vanishing PCAC quark mass is used. 

This constant, whose value in physical units is m<j> = 284(40) MeV, can be related 
to the quenched topological susceptibility xt through the Witten-Veneziano formula 
[39] 



f 



if^ 



R\2 



Xt , 



(37) 



which is valid in leading order of the l/N^ expansion. With xt = (193 it 9MeV)^ 
and our value for f^ we would obtain m$ = 426 MeV. 
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P=4-0 ^ 





0.01 0.02 0.03 0.04 0.05 

Figure 3: Pion masses squared in lattice units and the results of the 
PQChPT fit. 

5 Discussion 

This first Monte Carlo investigation of tlie hadron masses in QCD witli Nf = 1 dynam- 
ical quark flavour reveals the qualitative features of the low lying particle spectrum 
in this theory. The spatial extensions of our 12^ • 24 and 16^ • 32 lattices are about 
L ~ 2.2fm (see Table [2])0 This implies lattice spacings a ~ 0.19fm and a ~ 0.13 fm, 
respectively. The (bare) quark masses are reasonably small - in a range 10-30 MeV 
and 25-60 MeV on the 12^ • 24 and 16^ • 32 lattice, respectively. The updating algo- 
rithm we use (PHMC with stochastic correction [18]) works fine in this range making 
the extension of the Monte Carlo investigations towards larger volumes, smaller quark 
masses and smaller lattice spacings straightforward. In the present runs the fluctua- 
tion of the eigenvalues of the fermion matrix towards exceptionally small (or negative) 
values can be easily handled by reweighting the configurations during the evaluation 
of expectation values. In fact, except for the run with the smallest quark mass on the 
12^-24 lattice where the reweighting has a small effect, the reweighting is completely 
negligible or even unnecessary. 



^In order to have some relation to the scales in real QCD, we set the Sommer scale parameter by definition 
to be ro = 0.5 fm 
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Figure 4: Pion decay constants in lattice units and the results of the 
PQChPT fit. 



|3=3.8 

P=4-0 ^ 



The lightest hadron is the pseudoscalar meson bound state of a quark and an 
antiquark - the r^s-meson (see Tableland Figure [1]). The corresponding scalar bound 
state - the cxs-meson - is in our points by about a factor 1.5 heavier. Compared to 
the estimate in rn^^/mr^^ ~ Nd {Nc — 2) = 3 this result is too low but the situation 
could be better in the zero quark mass limit where the prediction of [8j applies to. 
The lightest baryon - the A^-baryon ~ is by a factor of about 3 heavier than the r/^- 
meson. The lightest glueball lies between the o"s-meson and the A^-baryon, but its 
mass could not be properly measured on the 16'^ • 32 lattice with our statistics. In 
general, the mass measurements have relatively large errors - between 3-10% ~ and 
no infinite volume and continuum limit extrapolations could be performed with our 
present data. We hope to return to these questions and to give more precise results in 
future publications. 

An interesting aspect of Nf = 1 QCD is the possibility of a partially quenched 
extension with valence quarks. In particular, adding two valence quarks, the model 
has similarities to QCD in nature with its three light (u, d and s) quark flavours. A 
theoretically interesting special case is if all three quarks, the dynamical one and the 
two valence ones, have exactly equal masses. In this case there is an exact SU(3) flavour 
symmetry. This can be exploited for the introduction of a quark mass by defining it 
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as the PC AC quark mass in the partiahy quenched theory. In this extended model 
there exist the usual light hadron states well known from real QCD: the pseudoscalar 
pseudo-Goldstone bosons (pions etc.), the nucleon etc. The results for the masses of 
the lightest states and the decay constant of the pseudoscalar bosons are collected in 
Table m and also shown in Figure [2j 

Since the physical volumes of the 12^ and 16^ lattices are to a good approximation 
equal, the comparison of the results at the two different lattice spacings gives a hint for 
the magnitude of the deviations from the continuum limit. As one can see in Figs. [1] 
andO the scaling between (3 = 3.8 and (3 = 4.0 is reasonably good - especially for the 
lightest states r]s and tt. However, for reliable continuum limit estimates more data at 
several lattice spacings are required. 

In the pseudoscalar sector of the partially quenched model one can apply partially 
quenched Chiral Perturbation Theory for fitting the mass and the decay constant. 
As Figs. E] and U] show, the NLO formulae give good fits but the number of degrees 
of freedom in the fits is small and therefore the uncertainty of the fit parameters is 
relatively large. 
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